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Abstract 



We analyze the critical behavior of two-dimensional iV-vector spin systems 
with noncollinear order within the five-loop renormalization-group approxi- 
mation. The structure of the RG flow is studied for different N leading to 
the conclusion that the chiral fixed point governing the critical behavior of 
physical systems with N = 2 and N = 3 does not coincide with that given by 
the 1/N expansion. We show that the stable chiral fixed point for N < N* , 
including N = 2 and N = 3, turns out to be a focus. We give a complete 
characterization of the critical behavior controlled by this fixed point, also 
evaluating the subleading crossover exponents. The spiral-like approach of 
the chiral fixed point is argued to give rise to unusual crossover and near- 
critical regimes that may imitate varying critical exponents seen in numerous 
physical and computer experiments. 
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I. INTRODUCTION 



The critical behavior of two-dimensional frustrated spin models with noncollinear or 
canted order has been the object of intensive theoretical and experimental studies, being of 
prime interest for those investigating layered magnetic systems with special structures and 
superconducting Josephson-j unction arrays in an external magnetic field. 

In physical magnets, frustration (leading to noncollinear order) may arise either because 
of the special geometry of the lattice or from the competition of different kinds of interactions. 
An example of the first kind is the two-dimensional triangular antiferromagnet, whereas for 
the latter case typical models are the fully frustrated XY model on a square lattice and 
models with different nearest-neighbor and next-nearest-neighbor interactions. In these 
systems the Hamiltonian is minimized by the noncollinear configurations with a 120° spin 
structure. As a consequence, at criticality, there is a breakdown of O(N) symmetry in the 
high-temperature phase to 0(iV — 2) symmetry in the ordered phase. 

Field-theoretical (FT) studies of systems with noncollinear order are based on the 
O(iV) xO(M) symmetric Landau- Ginzburg- Wilson (LGW) Hamiltonian [1-3] 

w = / \ E [( d M 2 + + (e 2 E [to • - M] } > (i) 

where tp a (1 < a < M) are M sets of iV-component vectors. We will consider the case M — 2, 
that, for vo > 0, describes frustrated systems with noncollinear ordering. Negative values of 
Vq correspond to simple ferromagnetic or antiferromagnetic ordering, and to magnets with 
sinusoidal spin structures [2]. 

The physically relevant cases are described by the Hamiltonian (1) with N = 2 (frus- 
trated XY model) and N = 3 (frustrated Heisenberg model). Despite the very intensive 
theoretical and experimental efforts to fully understand the critical behavior of these mod- 
els, the situation is still controversial. For the frustrated XY model the strongly debated 
point is either the single critical temperature T c exists at which both the SO (2) and the 
Z2 symmetries are simultaneously broken or there are two successive phase transitions at 
different critical temperatures. In the latter case, the order at which two phase transitions 
occur and the numerical values of the critical exponents are doubtful too. Obviously, study- 
ing the Hamiltonian (1), we cannot clarify the point just mentioned, since this Hamiltonian 
describes only the chiral phase transition. In more detail, these and other relevant issues 
are reviewed in Ref. [4], where one can find a complete list of references; we mention here 
only the very recent studies [5-8] (not quoted in Ref. [4]). 

The critical behavior of the frustrated Heisenberg antiferromagnet in two dimensions 
is clearer: it undergoes one phase transition mediated by Z 2 topological defects. Again, a 
complete list of references can be found in Ref. [4]. Here the debated point concerns whether 
the LGW Hamiltonian is able to keep the non trivial topological excitations present in this 
model, to the contrary of other approaches like the nonlinear a model (NLa). In Ref. [4] it 
is claimed that it may be the case and we will assume it in the following. 

The LGW Hamiltonian (1) has been extensively studied in the framework of e = 4 — d 
[2,9,10], and by means of the 1/N [2,10] expansion. The existence and the stability properties 
of the fixed points were found to depend on N and on the spatial dimensionality [2,3]. Within 
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the e expansion, for sufficiently large values of N the renormalization group equations possess 
four fixed points: the Gaussian fixed point (v — u — 0), the 0(2N) (v = 0) one and 
two anisotropic fixed points located in the region u, v > and usually named chiral and 
antichiral. The chiral fixed point is the only stable one. There is a critical dimensionality 
N c where the chiral and antichiral fixed points coalesce and disappear for N < N c . In the 
last case, under the absence of any stable fixed point the system undergoes a weak first-order 
phase transition, since the associated RG flows run away from the region of stability of the 
fourth-order form in the free energy expansion. The three-loop estimate of iV c is [9] 

iV c = 21.8-23.4e + 7.1e 2 + 0(e 3 ), (2) 

that, after an appropriate resummation, results in iV c > 3 in three dimensions. This inequal- 
ity leads to the conclusion that for the physical models with N = 2, 3 the three-dimensional 
chiral transition is first order, as corroborated also by some other RG studies [11,12]. To 
the contrary, both for N = 2 and N = 3 the highest-order six-loop calculations in three 
dimensions [13] reveal a strong evidence for a stable chiral fixed point that, however, is not 
related to its counterpart found within the e expansion. The four-loop two-dimensional anal- 
ysis shows an equivalent topology of the fixed points location in the renormalized coupling 
constants plane [4]. Furthermore, in a recent work [14] we claimed that both in two and 
three dimensions the controversial situation may reflect the quite unusual mode of critical 
behavior of the iV-vector chiral model under the physical values of N. It was shown that 
this critical behavior is governed by a stable fixed point which is a focus, attracting RG 
trajectories in a spiral-like manner. Approaching the fixed point in a nonmonotonic way, 
looking somewhat irregular, may result in a large variety of the crossover and near-critical 
regimes. 

In this paper, we report the perturbative five-loop two-dimensional RG expansions and 
carry out their analysis. Some of the results for N = 2 and N = 3 was already anticipated 
in [14]. We point out the idea that the critical behavior for large N is the one predicted by 
the e expansion. With decreasing N, the stable chiral fixed point becomes focuslike at some 
marginal value of N, N* that differs, in principle, from N c where both chiral and antichiral 
fixed points disappear. With the precision of our calculations, we are not able to fully clarify 
if there is a region of values of N, where the chiral transition is first order. It is worthy 
to note that the parallel work on the three-dimensional models [15] reveals a full analogy 
between chiral critical behaviors in two and three dimensions. 

The paper is organized as follows. In Sec. II we derive the perturbative series for the 
renormalizat ion-group functions up to five-loop order. The resummations methods and the 
results of the analysis are presented in Sec. III. In Sec. IV we draw our conclusions. 

II. PERTURBATIVE EXPANSIONS IN TWO DIMENSIONS 

A. Renormalization of the theory 

The fixed-dimension field-theoretical approach [16] represents an effective procedure in 
the study of the critical properties of systems that undergo second-order phase transitions 
[17]. In the case under consideration, the expansion is performed in two quartic coupling 
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constants of the Hamiltonian (1). The theory is renormalized by a set of zero-momentum 
conditions for the one-particle irreducible two-point correlation function, four-point correla- 
tion function, and two-point correlation function with an insertion of the operator |0 2 : 

,(2) 



rg w (p) = M,,^ 1 m 2 + P 2 + o(p 



(3) 



ritL(o) = V™ 



u 



— (5 a b5 c d + S ac Sbd + Sadhc) + V C a i,bj,ck,dl 



(4) 



rii ,2) (o) = ^z t - 1 , 



(5) 



where C a i^j, c k,di is an appropriate combinatorial factor [4,13]. 

The perturbative knowledge of the functions r*- 2 -*, r'- 4 *', and r'- 1 ' 2 -' allows one to relate the 
renormalized physical parameters (u,v,m) to the bare ones (uo,vo,r ). 

The fixed points of the model are defined by the common zeros of the (3 functions, 



fJ u {u,v) = m — 



u ,v 



P v (u,v) = m^ 



(6) 



u ,v 



The stability properties of these points are determined by the eigenvalues Ui of the matrix: 



ft 



( dp u {u,v) dp u (u,v) \ 



du dv 
d(3 v (u,v) d(3 v (u,v) 



V du 



(7) 



dv 



A fixed point is stable if both the eigenvalues have a positive real part. If the eigenvalues 
possess nonvanishing imaginary parts, the fixed point is called a focus and the corresponding 
RG trajectories are spirals. The eigenvalues uj,i are connected to the leading scaling correc- 
tions, which go as ~ |t| Al = \t\ VWi , with few exceptions as the two-dimensional Ising 
model (see for a discussion Ref. [18]). 

The values of the critical exponents r], v and 7 are related to the RG functions and 
rj t evaluated at the stable fixed point: 



77 = rj^u ,v ), 



(8) 



where 



u=[2- V<p (u*,v*)+Vt(u*,v*)]- 1 , 



r)4>{u,v) = 



d\nZ, 



dlnm 



d\nZ^ d\nZ, 

Pu rs. \ Pv 



Uo,VQ 



du 



dv 



(9) 



(10) 



r} t (u,v) = 



d la Z f 



dlnm 



d\nZ t d\nZ t 

Pu 7\ \~ Pv' 



UQ,V0 



du 



dv 



(11) 
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B. Five-loop series 



In this section we extend the four-loop perturbative series [4] for the RG functions fi u , 
/3 V , r)<f„ and r) t up to five- loop order using the numerical values of two-dimensional five-loop 
integrals evaluated in Ref. [19]. 

In order to obtain finite fixed point values in the limit of infinite components of the order 
parameter (N — > oo) we use the rescaled couplings 

8tt „ 8tt n 

u = — R 2 n u, v = — R 2N v, (12) 

where R N — 9/(8 + N). 

The resulting RG series are [20]: 

* = +* 2 + <^r - wrm* + (13) 

- 1 6-iV n 6 / n . . 

* = " 2 (IT^f + ^TN) UV + ^.^ ' ( } 

ru ~ L83417(1 + 7V) ^ + 1-83417 (1-JV)" _ 1.37563 (1 - AQ w 

^" (8 + 2iV)2 M+ (8 + 2iV)2 MV (8 + 2iV)2 + ^ >3 ij ' [ j 

_2(1 + AQ_ (l-iV)_ 13.5025(1 + AO _ 2 13.5025 (1 - iV) 
^ ~ " (4 + iV) U ~ JA + N) V + (8 + 2Nf U + (8 + 27V) 2 UV 
_ 10.1269(1 -AQ w 

(8 + 2iV) 2 + J+ t>3 " ' ( j 

where 

- 3 i 3 i 

flu = R 2 N fiui fiv = R 2 N0v (17) 

167T 167T 

Only the five-loop coefficients b[j\ e\f and are reported in Table I, since up to the 
fourth order they are reported in Ref. [4]. Note that, due to the last rescaling, the element 
of the Q matrix (7) are two times the derivative of j3 with respect to u and v. 

We verify the exactness of our series by checking them in particular limits (N = 1, v = 0, 
N = M = 2, N -> oo) as in Ref. [4]. 



III. RESUMMATION AND ANALYSIS OF THE FIVE-LOOP SERIES 

The field theoretical perturbative expansions are known to be divergent asymptotic se- 
ries. Quantitative informations may be extracted from these series exploiting the properties 
of Borel summability for 4 theories in two and three dimensions and resumming them by 
a Borel transformation combined with a method for the analytical extension of the Borel 
transform. This last procedure can be realized by using Pade approximants or via a con- 
formal mapping, which maps the domain of analyticity of the Borel transform (cut at the 
instanton singularity) Ub onto a circle [21]. 
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TABLE I. Five-loop coefficients e$ and ef) of RG functions. 



i i R~ J b 



3lW 



+2.01674 + 2.10643 N + 0.716897 N 2 + 0.0887555 N 3 +0.00241689 JV 4 - 2.45381 * 10 _y N b 

5.1 +2.92727 - 1.33251 N - 1.34924 N 2 - 0.238028 N 3 -0.00748703 N 4 - 0.000013119 N 5 

4.2 -4.43837 + 2.03796 N + 2.0376 N 2 + 0.352602 N 3 +0.0101649 N 4 + 0.0000459364 AT 5 

3.3 +3.19454 - 1.49238 N - 1.45358 N 2 - 0.242665 N 3 -0.0058521 N 4 - 0.0000594059 TV 5 

2.4 -1.15807 + 0.561265 JV + 0.515855 N 2 + 0.0797969 N 3 +0.00112123 N 4 + 0.0000298703 N 5 

1.5 +0.222479 - 0.113395 N - 0.095816 N 2 - 0.0132978 N 3 +0.0000283847 N 4 + 1.2973 * 10~ 6 N 5 

0,6 -0.0151778 + 0.00751941 N + 0.00660415 N 2 + 0.00104974 N 3 +8.90815 * 10~ 6 N 4 - 4.38647 * 10~ 6 N 5 



5X1 +6.24588 + 3.59449 N + 0.620789 N 2 + 0.0264763 N 3 -0.000525164 JV 4 - 0.0000277103 N° 

4.2 -8.7215 - 4.66521 N - 0.686267 N 2 - 0.0125189 A^ 3 +0.00131132 N 4 + 0.0000937132 TV 5 

3.3 +6.32764 + 3.38603 N + 0.514699 N 2 + 0.0164667 N 3 -0.00034289 N 4 - 0.000120039 N 5 

2.4 -2.44897 - 1.35035 N - 0.227251 N 2 - 0.0119867 N 3 -0.0000336496 N 4 + 0.0000731821 N 5 

1.5 +0.48559 + 0.257994 N + 0.0396213 N 2 + 0.000895835 N 3 -0.000205336 N 4 - 0.0000218222 N 5 

0,6 -0.036149 - 0.0191764 N - 0.00232213 N 2 + 0.000377994 N 3 +0.0000915813 N 4 + 2.59344 * 10~ 6 TV 5 



h3 



r-«-J 
2N 



Li_ 



5,0 
4,1 
3,2 
2,3 
1,4 
0,5 



-0.00722954 - 
-0.0180738 + 
+0.0260483 - 
-0.0161648 + 
+0.00411472 - 
-0.000397627 



- 0.011019 N - 0.00391279 N 2 - 0.000142581 N 3 -0.0000192165 N 4 
0.00860027 N + 0.00916515 N 2 + 0.00026037 N 3 +0.0000480413 N 4 
0.0129479 N - 0.0129842 N 2 - 0.0000681925 N 3 -0.0000480413 N 4 
0.00849874 N + 0.00785318 N 2 - 0.000211121 N 3 +0.0000240207 N 4 

- 0.00235806 N - 0.00191523 N 2 + 0.000164576 N 3 -6.00516 * 10~ 6 N 4 

+ 0.000247705 N + 0.000181639 N 2 - 0.0000311162 N 3 - 6.00516 * 10" 7 N 4 



l 'J n 2N e ij 

5j0l -0.285954 - 0.463731 N -0. 195487 N 2 - 0.0176704 N 3 +0.000039754^* 

4.1 -0.714885 + 0.270442 N + 0.400167 N 2 + 0.0443747 N 3 -0.000099385 N 4 

3.2 +0.954626 - 0.384666 N - 0.520913 N 2 - 0.0491461 N 3 +0.0000993852 N 4 

2.3 -0.572386 + 0.253092 N + 0.29846 N 2 + 0.0208835 N 3 -0.0000496928 N 4 

1.4 +0.149502 - 0.0775891 N - 0.0701865 N 2 - 0.00179722 N 3 +0.0000708398 N 4 
0,5 -0.0157319 + 0.0101789 N + 0.00583765 N 2 - 0.00022754 N 3 -0.0000571742 N 4 



Let us consider a perturbative series in u and v 

l l—k 

R(u,v)= Y,Y, R hkU h v\ (18) 

fc=0 h=0 

which we want to resum (i.e., one of the 77 or j3 functions). In the Pade-Borel method we 
use the trick of the resolvent series considering the Borel transformed function, 

00 \n n 

R(X, u, v,b)= w n , un E Ri,n-iu l v n ~ l . (19) 
n =o L\n-rO-\- 1) l=Q 

Calling P^(X,u,v,b) the [N/M] Pade approximant of the series R(\,u,v,b) in the variable 
A, one obtains the estimate of the desired quantity as 

Pj*(b,u,v)= / e-H b P^(t,u,v,b)dt. (20) 
Jo 

In this way we produce several approximants of the function R(u,v) with varying the pa- 
rameters b, N, and M. 

In the conformal mapping method, the idea is to exploit, in the course of the resumma- 
tion, the knowledge of the large order behavior of the series R(u, v) at fixed z = v/u [22,13,4]. 
This large-order behavior is related to the singularity of the Borel transform closest to the 
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origin [23]. For the Hamiltonian (1), the value of the Borel transform singularity closest to 
the origin uj, has been computed in Ref [4], obtaining: 

-L = - a R 2N for < z < 4, (21) 

u b 

— = — a R 2 n ( 1 z ) for z < 0, z > 4, 



u b V 2 

where a = 0.238659217. . .. For z > 2 there is a singularity on the real positive axis which, 
however, is not the closest one to the origin for z < 4. Thus, for z > 2 the series are not 
Borel summable. We resum the RG functions also for the values of z for which the series 
are not Borel summable since the method should provide a reasonable estimate if z < 4, 
because we take into account the leading large-order behavior. 
With the knowledge of Ttb, one can perform the mapping 



y{u; z) = v = , (22) 

yjl - U/U b {z) + 1 

in order to extend R(u, uz) to all positive values of u. Consequently one obtains the set of 
approximants 

E(R) p (a,b;u,v) = J2B k (a,b-,v/u) x f°° dtt^ K , (23) 



fc=0 



[l-y(ut;v/u)Y' 



where the coefficient B\. is determined by the condition that the expansion of E(R) p (a, b; u, v) 
in powers of u and v gives R(u, v) to order p. 

An important issue in the fixed dimension approach to critical phenomena concerns 
the analytic properties of the f3 functions. As shown in Ref. [18] for the O(N) model, 
the presence of confluent singularities in the zero of the perturbative j3 function causes a 
slow convergence of the resummation of the perturbative series to the correct fixed-point 
value. The apparent stability of the results when analyzing a finite number of terms of the 
perturbative expansion may not provide a reliable indication of the uncertainty of the overall 
estimates. The O(N) two-dimensional field-theory valuations of physical quantities [21,19] 
are less accurate than the three-dimensional counterparts, due to the bigger critical value 
of the quartic coupling constant and stronger nonanalyticities at the fixed point [18,24,25]. 
In Ref. [18] it is explicitly shown that the nonanalytic terms cause a large deviation in the 
estimate of the right correction to the scaling, i.e., the exponent uj. At the same time, the 
perturbative results for the fixed-point location turn out to be rather good approximations 
for the exact ones. We think that a similar situation holds for frustrated models. 



A. Fixed points, stability, and critical exponents 

In order to calculate the fixed points of the theory, we resum the perturbative expression 
for each f3 function in the whole plane (u, v) following the procedure explained above. With 
the conformal mapping method, we choose the approximants that stabilize the series for 



7 




FIG. 1. Zeros of (5 functions for N = 2 in the (it, v) plane. Pluses (+) and crosses (x) correspond 
to zeroes of /3v(u,v) and /3u(u,v), respectively. 




FIG. 2. Zeros of f3 functions for N = 3 'm the (u, v) plane. Pluses (+) and crosses (x) correspond 
to zeroes of flv(u,v) and (3u(u,v), respectively. 

small values of the coupling constants with varying the considered perturbative order, i.e., 
a = 0,1,2 and b = 5, 7, . . . , 15. In this manner we have 18 approximants for each (3 
function. To obtain reasonable values and error bars for the fixed-point coordinates, we 
divide the domain 0<w<4, 0<17<6in 40 2 rectangles, and we mark all the sites in 
which at least two approximants for (3^ and (3$ vanish. This procedure is applied for the 
four-loop and five-loop series and the results for the zeros of the (3 functions for iV = 2 
and N = 3 are displayed in Figs. 1 and 2, respectively. We remind the reader that the 
four- loop results were previously reported in Ref. [4], where similar figures were shown also 
for the three-loop case. It was also shown that no upper branch of zeros of (3^ appears in 
the three-loop approximation for all considered values of N. 

As one can see from Figs. 1 and 2, for these physical systems the resummed RG expan- 
sions yield similar pictures for the zeros of the (3 functions. The existence of three fixed 
points - the gaussian, the 0(2N), and the chiral - is clear in both cases, while the presence 
of an unstable antichiral fixed point looks doubtful. In fact, if this last point exists, it should 
be located in the domain where the resummation procedure fails to give reliable quantitative 
results (z > 4). 

The location of the chiral fixed point for N = 2 is 

[u*,v*] = [2.3(2), 5.0(5)] (4-loopCM), (24) 

[2.25(25), 4.80(45)] (5-loop CM), (25) 

[2.35(30), 4.35(25)] (Pade), (26) 

and for N = 3 
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[u*,v*]= [2.3(3), 3.9(5)] (4-loop CM), 

[2.2(2), 3.9(3)] (5-loopCM), 
[2.25(20), 3.60(15)] (Pade). 



(27) 
(28) 
(29) 



The stability properties of the chiral fixed point depend on the signs of the real parts 
of the eigenvalues u>i of the Q matrix. We take the above 18 approximants for each f3 
function and compute the numerical derivatives of each couple of approximants of the two 
f3 functions at their common zero, obtaining 324 possible combinations. To have reliable 
numerical estimates for u>i, we limit ourselves by the approximants that yield the chiral fixed 
point coordinates compatible with their final, properly weighted values. 

Both for N = 2 and iV = 3, we find that the chiral fixed point is a stable focus since 
it possesses stability eigenvalues with nonvanishing imaginary parts and positive real parts. 
At the five-loop (four-loop) approximations for N = 2 the complex eigenvalues are produced 
by 87% (87%) of the working approximants within the conformal mapping method, while 
for N = 3 by 89% (74%) of the approximants employed. When estimating the averaged 
eigenvalues of the Q matrix at the chiral fixed point, we disregard the approximants leading 
to purely real eigenvalues. For N = 2 the final results of our calculations are: 



and for N = 3 



u± = 1.50(25) ± i 1.00(45), (4-loop CM) 
u± = 2.05(35) ± i 0.80(55), (5-loop CM) 



u± = 1.30(25) ± i 0.50(35), (4-loop CM) 
u± = 1.55(25) ± i 0.55(35), (5-loop CM) 



(30) 
(31) 



(32) 
(33) 



The presence of imaginary parts comparable in magnitude with the real ones leads to 
the conclusion that the critical behaviors of frustrated two-dimensional Heisenberg and XY 



TABLE II. Location of the chiral fixed point with its stability eigenvalues, critical exponent 77, 
and crossover exponent u>o(2N) a t the 0(2N) fixed point. 



N = 2 


Pade 


CM. (4 loop) 


CM. (5 loop) 


— * 

V 

u± 

7] 

w O(4) 


2.35(30) 
4.35(25) 

2.33(65) ±i0.62(48) 
0.3(1) ' 


2.3(2) 
5.0(5) 

1.50(25) ±il.00(45) 

0.29(5) 

-0.36(4) 


2.25(25) 
4.80(45) 

2.05(35) ±i0.80(55) 

0.28(8) 

-0.35(16) 


N = 3 


Pade 


C.M.(4 loop) 


CM. (5 loop) 


u* 

— * 

V 

u± 

7] 

w O(6) 


2.25(20) 
3.60(15) 

1.77(43) ±»0.47(18) 
0.23(7) 


2.3(3) 
3.9(5) 

1.30(25) ±i0.50(35) 

0.24(6) 

-0.7(2) 


2.2(2) 
3.9(3) 

1.55(25) ±i0.55(35) 

0.23(5) 

-0.69(9) 
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u 



1 2 3 4 1 2 3 4 



FIG. 3. Zeros of (3 functions for N = 32 in the (u,v) plane. Pluses (+) and crosses (x) 
correspond to zeroes of (3v(u,v) and (3u(u,v), respectively. 




FIG. 4. Zeros of (3 functions for N = 16 in the (u,v) plane. Pluses (+) and crosses (x) 
correspond to zeroes of (3\j(u,v) and (3u(u,v), respectively. 

systems are driven by focus fixed points. To strengthen this new and important issue, we 
repeat the analysis of the stability properties of the chiral fixed point by a Pade study, which 
confirms the previous results, although the statistics in this case is much less significant due 
to a large number of defective Pade approximants. All the results obtained are summarized 
in Table II. Note that in the past such a behavior was found in three dimensions within the 
three-loop RG approximation but only for unphysical values of iV (N = 5, 6, 7) [27]. 

In order to fully characterize the two physical systems discussed, we calculate the critical 
exponents as well. For the exponent 77 we find values that agree with the results of the 
four-loop analysis [4]. At the same time, we can not give a reasonable estimate for the 
critical exponent v because of the strong oscillations observed with varying the working 
approximants. 

Finally, we report in Table II the five- loop estimate of the crossover exponent lvo(2N), 
finding that the Heisenberg fixed point is unstable for both N = 2, 3, in agreement with the 
previous four- loop calculation [4]. 

We apply the same kind of analysis also to the models with N > A. All these systems are 
free of topological defects [4] and are well described by appropriate NLa models [26] which 
provide the existence of a stable chiral fixed point with critical exponents and corrections 
to the scaling equal to the ones of a system with an infinite number of components of the 
order parameter. 

The zeros of the two (3s for N = 32, 16, 8, 6, 4 are shown in Fig. 3, 4, 5, 6, 7. From these 
figures it is clear that for large enough values of N four fixed points exist: the gaussian, the 
0(2N), the chiral, and the antichiral, with a topological structure similar to the one found 
in e expansion for N > N c . The location of the chiral fixed point for various N is reported 
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FIG. 5. Zeros of (3 functions for N = 8 in the (it, v) plane. Pluses (+) and crosses (x) correspond 
to zeroes of (3\jiu, v) and /3u(u, v) respectively. 



v 41 v 51 




FIG. 6. Zeros of f3 functions for N = 6'm the (u, v) plane. Pluses (+) and crosses (x) correspond 
to zeroes of /3y(u, v) and Pu(u, v) respectively. 




FIG. 7. Zeros of (3 functions for N = 4 in the (u, v) plane. Pluses (+) and crosses (x) correspond 
to zeroes of Pv(u, v) and 0u(u, v) respectively. 

in Tab. III. 

From the NLa model it is expected lo\ = uj 2 = 2 [26]. For this reason iV* should 
be less or at least equal to 4. From our analysis we find that the chiral fixed point is 
stable for all the considered N. However, the estimates for uo,i are slighty different from 
2, also for large N (e.g., for iV = 16 we find ^ = 1.59(5) and cu 2 = 2.01(11)). Such a 
disagreement was already observed in the four-loop approximation, and it may be attributed 
to the nonanalyticities discussed in the previous subsection. We observe that for iV > 8 all 
the stability eigenvalues are purely real, but for N = 6 about the 50% of the five-loop 
approximants possess a nonvanishing imaginary part, suggesting iV* ~ 6. This result seems 
to disagree with the NLa prediction. Again this might be due to nonanalyticities. 

To further approve the above conclusions, we apply, along with the conformal-mapping 
based, the resummation by means of the Pade-Borel technique. In the region where the 
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TABLE III. Location of the stable chiral fixed point for N > 4. 



N 


C.M.(4 loop) 


C.M.(5 loop) 


Pade 


4 




[2.2(1),3.37(23)] 


[2.22(16),3.19(10)] 


6 


[2.26(7),2.85(15)] 


[2.21(6),2.80(8)j 


[2.25(10),2.79(4)] 


8 


[2.28(9),2.66(8)] 


[2.22(7) ,2.61(6)] 


[2.23(6),2.61(3)] 


16 


[2.22(4),2.37(3)j 


[2.17(5),2.33(3)j 


[2.18(2),2.33(3)j 


32 


[2.128(16),2.19(1)] 


[2.09(2),2.16(2)j 


[2.100(15),2.165(10)] 



fixed point is expected, many Pade approximants turn out to be defective, i.e., possessing 
dangerous poles. We find that the approximants [4/1], [3/1], and [2/2] are nondefective in 
the majority of the cases. Considering these three Pade approximants with b = 0, 1, 2 (being 
the more stable under variation of the number of loops), we finally have nine approximants 
for each (3 function. The location of the chiral fixed point is reported in Table III, with all 
the defective cases discarded. 

We also compute the exponents rj and i] t with both the resummation techniques. We 
consider a lot of approximants and finally we choose the more stable ones with varying 
the number of loops. The found values are close to those of the O(N) model for N > 3, 
but differ from the exact ones coming from NLa model i] = and r) t = —2 (e.g., we find 
rjt = —1.95(20) for N = 16 and rj = 0.16(2) for N = 6). We think these discrepancies 
are due to nonanalyticities, as in the case of O(N) models. To conclude the analysis of the 
fixed points, we study the stability of the Heisenberg fixed point by calculating the crossover 
exponent u>o(2N)- The five- loop results show that this fixed point is unstable, in agreement 
with the conclusion drawn earlier from the four- loop series [4]. 




FIG. 8. RG flow for the two-dimensional chiral model with N = 2. The approximants used are 
given by a = 1 and b = 5 for j3u and a = and b = 8 for fa. In this case, the chiral fixed point is 
at (2.427,5.045). 
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0.5 1 1.5 2 



FIG. 9. Crossover of the effective exponent rj e ff for N = 2. 

B. Renormalization-group flow and crossover 

In this subsection we demonstrate the structure of the RG flow for physical values of N. 
As it was already shown, a stable focus governs the critical behavior of these systems. In 
Fig. 8 (N = 2) one can see the spiral-like approach of the stable fixed point which is peculiar 
for the fixed points with stability eigenvalues possessing nonvanishing imaginary parts. The 
shaded area is the region where the resummation procedure is expected to fail (i.e., z > 4), 
since we do not take into account the singularity of the Borel transform closest to the origin 
that is on the real positive axis (cf. Eq (21)). Consequently, the RG trajectories given by 
the resummed perturbative series are quantitatively correct only within the unshaded areas. 
We report, nevertheless, the RG flows in other parts of the coupling constant plane in order 
to present a complete qualitative picture. Note, that the domain usually referred to as the 
sector of the first-order transitions, where the positiveness of the quartic form in the free- 
energy expansion breaks down, is given by the inequality z > 2. In fact, however, because of 
the presence of the higher-order terms in this expansion keeping the system globally stable 
at any temperature, the true domains of the first-order transitions should be more narrow. 

From Fig. 8 we see that for some initial values of (u,v) the RG trajectories have the 
coordinate v that grows very fast at the beginning and seems to reach the region of the first- 
order phase transitions, but just before arriving there these trajectories drastically curve 
moving toward the stable chiral fixed point. These bizarre trajectories can imitate both the 
discontinuous phase transitions and strongly scattered effective critical exponents observed 
in numerous physical and computer experiments. To substantiate the former scenario, we 
calculate the numerical values which the effective exponent r] e ff takes along the RG trajec- 
tories running from the Gaussian fixed point to the chiral one. Corresponding curves are 
shown in Fig. 9. As is seen from this figure, the exponent oscillates within a large range, 
even near the stable fixed point. 

We demonstrate also the RG flows for iV = 3 (Fig. 10) and for iV = 8 > iV* (Fig. 
11). In the former (physical case) the stable fixed point is a focus and phase trajectories 
approach it in a spiral-like manner. In the latter (unphysical case) the stable fixed point 
turns out to be a node and the trajectories are certainly not spirals. Also for N = 3 we 
estimate the effective exponent r) e ff along the RG trajectories running from the Gaussian 
fixed point to the chiral one. The corresponding curves (cf. Fig. 12) show big oscillations 
and sometimes reach negative values. We note that these negative values are obtained only 
in the region where the resummation is not sure, i.e., z > 4. Nevertheless, this strange 
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FIG. 10. RG flow for the two-dimensional chiral model with N = 3. Both (3 functions are given 
by approximants given by a = 2 and 6 = 9, the chiral fixed point is located at (1.702, 2.858). 




0.5 1 0(16)1.5 2 2.5 3 

FIG. 11. RG flow for the two-dimensional model with N = 8. The approximants used are given 
by a = and b = 7 for f3u and a = 1 and b = 5 for f3y. In this case, the chiral fixed point is at 
(2.226,2.571) and the antichiral fixed point is at (0.891,3.133). 



crossover, if reproduced from the three-dimensional series [15], might explain some doubtful 
experimental results, where a negative rj was found. 

To conclude, it is worthy to note that earlier the focuslike stable fixed points were found 
on the RG flow diagrams of the model describing critical behavior of liquid crystals [28] and 
of the 0(n)-symmetric systems undergoing first-order phase transition close to the tricritical 
point [29] . In those cases, however, the independent coupling constants had different scaling 
dimensionality and played essentially different roles in critical thermodynamics. Moreover, 
recently complex subleading exponents were found from the high-temperature expansion of 
Dyson's hierarchical Ising model [30]. 



IV. CONCLUSIONS 

We studied the critical thermodynamics of the two-dimensional iV-vector chiral model 
in the five-loop RG approximation. Using the advanced resummation technique based upon 
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FIG. 12. Crossover of the effective exponent rj e ff for N = 3. 



Borel transformation combined with conformal mapping and use of Pade approximants, it 
was shown that the fixed point governing the chiral critical behavior is a stable node only for 
large enough N, N > N*. For smaller N a chiral fixed point turns out be a stable focus and 
the RG trajectories approach this point in a spiral-like manner. Of particular importance is 
the fact that a focus-driven critical behavior was found for systems with physical values of N, 
N = 2, 3. Evaluating the critical exponent 77 along some RG trajectories, we demonstrate 
that spiral-like approach of the chiral fixed point results in unusual crossover and near- 
critical regimes that may imitate varying critical exponents seen in physical and computer 
experiments. 
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